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Abstract 

Recently, the DAMA/LIBRA collaboration has repeated and reinforced their claim to have detected 
an annual modulation in their signal rate, and have interpreted this observation as evidence for dark- 
matter particles at the 8.2<r confidence level. Furthermore, it has also been noted that the effects of 
channeling may enable a WIMP that scatters elastically via spin-independent interactions from nuclei 
to produce the signal observed by DAMA/LIBRA without exceeding the limits placed by CDMS, 
XENON, CRESST, CoGeNT and other direct-detection experiments. To accommodate this signal, 
however, the mass of the responsible dark-matter particle must be relatively light, moM ^ 10 GeV. 
Such dark-matter particles will become captured by and annihilate in the Sun at very high rates, 
leading to a potentially large flux of GeV-scale neutrinos. We calculate the neutrino spectrum resulting 
from WIMP annihilations in the Sun and show that existing limits from Super-Kamiokande can be 
used to close a significant portion of the DAMA region, especially if the dark-matter particles produce 
tau leptons or neutrinos in a sizable fraction of their annihilations. We also determine the spin- 
dependent WIMP-nuclei elastic-scattering parameter space consistent with DAMA. The constraints 
from Super-Kamiokande on the spin-dependent scenario are even more severe — they exclude any self- 
annihilating WIMP in the DAMA region that annihilates 1% of the time or more to any combination 
of neutrinos, tau leptons, or charm or bottom quarks. 
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I. INTRODUCTION 



Recently, the DAMA collaboration has provided further evidence for an annual modulation in 
the rate of nuclear- recoil events observed in their experiment [l|. Such a signal arises naturally 
from postulating weakly interacting massive particles (WIMPs) in the Galactic halo that scatter 
from target nuclei in detectors. The annual modulation of the interaction rate results from 
the variation in the relative velocity of the Earth with respect to the Galactic dark-matter 
halo as the Earth orbits the Sun. This changes the flux of dark-matter particles and their 
velocity distribution, with expected extrema occurring at June 2 and December 2. The DAMA 
experiment observes a maximum rate at low nuclear-recoil energies on May 24, plus or minus 
8 days, and have accumulated enough data to put the significance of the observed modulation 
at approximately 8<r. Both the phase and amplitude of the signal are highly suggestive of 
WIMP interactions. The collaboration has not been able to identify any other systematic 
effects capable of producing this signal, and have claimed that the annual modulation is a 
discovery of dark matter. This claim has been controversial, partly because a number of other 
experiments appear to be in direct contradiction. 

Several studies have attempted to reconcile the DAMA modulation signal with the null 
results of other direct-detection experiments 0, H, 0, d, @, 0], assuming a spin-independent 
WIMP-nucleus interaction. A feature common to the analyses of Refs. 0, 0, 0, @, H, 0] is that 
an elastically scattering WIMP with a mass in the several GeV range can satisfy the results 
of DAMA. The consistency of this region with the null results of CDMS [10], CRESST El 



CoGeNT [12J , and XENON [13j depends on the details of how the DAMA recoil energy spectrum 
is fit. If the data is divided into the two bins 2 — 6 keVee and 6 — 14 keVee, where keVee is the 
electron-equivalent recoil energy, then DAMA and the null experiments can be simultaneously 
accommodated by 3— 8 GeV WIMPs. If the DAMA modulation data is binned more finely, then 
the modulation in the 2 — 2.5 keVee bin right at the DAMA threshold is difficult to reconcile 
with other direct-detection constraints. We refer the reader to the studies in Refs. H, 0] for 
further details. The allowed parameter region depends crucially on the occurrence of channeling 



in the Nal crystals of the DAMA apparatus, an effect noted in Ref. [14| and studied by the 



DAMA collaboration 15 



In this paper, we consider the constraints to the new DAMA dark-matter parameter space 



that come from null searches for energetic neutrinos from the Sun J16j. If WIMPs scatter 
from nuclei in DAMA, then they will also scatter from nuclei in the Sun, be captured and 
thus annihilate therein, and thus produce energetic neutrinos that can be sought in experi- 
ments such as IceCube [T^, AMANDA [3], Baksan 0, MACRO Q, ANTARES [Hj, and 
Super-Kamiokande 22|. The only caveat is that the energetic- neutrino spectra will depend 
on the WIMP-annihiation products. Still, the range of neutrino spectra is bracketed, and so 



model-independent bounds can be sometimes obtained [23]. In particular, the Sun is composed 



primarily of protons (nuclei with spin), and so null neutrino searches should be especially con- 
straining for an explanation of DAMA in terms of a WIMP with spin- dependent interactions, 
as shown in Ref. 24j for the old spin-dependent DAMA parameter space. 

Here, we determine the regions of the mass-cross-section parameter space implied by the 
new DAMA results for a spin-dependent WIMP, and we determine parameter space (for both 
spin-dependent and spin-independent WIMPs) eliminated by null neutrino searches. In order 
to extend the analysis of Ref. [24] to the lower WIMP masses implied by channeling, we calcu- 
late the full neutrino energy spectra produced by decays of tau leptons and charm and bottom 
quarks, and we consider prompt annihilation to neutrinos. We also include the effects of neu- 
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trino mixing and the effects of WIMP evaporation from the Sun. The muon energy thresholds 
of ANTARES and IceCube are 10 GeV or higher and are therefore not able to observe the 
neutrinos produced by a several-GeV WIMP. We therefore consider bounds arising from the 
Super-Kamiokande experiment, which can identify muons with energies as low as 1.6 GeV. The 
limits imposed by Baksan and MACRO are similar to those coming form Super-Kamiokande. 
We study the DAMA allowed region arising from both the two-bin study of Ref. [E| and the 
parameter space arising from the full spectral analysis both with and without the 2 — 2.5 keVee 
bin 

Our conclusions are that that upper limits to the flux of energetic neutrinos from the Sun 
severely constrain the DAMA spin-independent parameter space if WIMPs annihilate directly 
to neutrinos or to tau leptons. The constraints from Super-Kamiokande on the spin-dependent 
scenario are even more severe; they exclude any spin-dependent WIMP in the DAMA region 
that annihilates 1% of the time or more to any combination of neutrinos, tau leptons, or charm 
or bottom quarks. 

Our paper is organized as follows. We review in Section [Til the formalism of WIMP capture 
and annihilation in the Sun, including the effect of evaporation, which is important for dark 
matter near the 3-GeV lower edge of the DAMA allowed region. In Section [TTTI we discuss the 
detection of the neutrinos from WIMP annihilation using upward-going muons in the Super- 
Kamiokande detector and study the constraints imposed by this process on the DAMA allowed 
region. We also present in Sec. IHII the spin-dependent parameter space (including channeling) 
implied by DAMA. In Sec. IIV} we comment on the energetic-neutrino constraint as applied to 
the case of light neutralino dark matter. Finally, we summarize our conclusions in Section |Vj 



II. WIMP CAPTURE AND ANNIHILATION IN THE SUN 



We briefly review here the basic formulae describing WIMP capture and annihilation in the 
Sun. A generic species of dark-matter particle present in the solar system will scatter elastically 
with and become captured in the Sun at a rate given by [25] 
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where pdm is the local dark-matter density, v is the local root-mean-square velocity of halo 
dark-matter particles, and is the dark-matter mass. <th and an e are the elastic scattering 
cross sections of the WIMP with hydrogen and helium nuclei, respectively. The factor of 1.1 
reflects the solar abundance of helium relative to hydrogen and well as dynamical factors and 
form factor suppression. The quantity S is a kinetic suppression factor given by 
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For WIMPs much heavier than hydrogen and helium nuclei, this suppression can be consider- 
able. For WIMPs in the 1-10 GeV range, however, S rs 0.9 — 1.0. 

If the capture rate and annihilation cross sections are sufficiently large, equilibrium will be 
reached between these processes. The differential equation governing the number of WIMPs in 
the Sun, denoted here as N, is 

N = C & — A G N 2 - E G N, (4) 

where C & is the capture rate, A Q is the annihilation cross section times the relative WIMP 
velocity per unit volume, and E & is the inverse time for a WIMP to exit the Sun via evaporation. 
A & can be approximated by 

a- = M, (5) 

VeE 

where V e s is the effective volume of the core of the Sun determined roughly by matching the 
core temperature with t he g ravitational potential energy of a single WIMP at the core radius. 
This was found in Refs. [2a, |27j to be 

V eS = 5.7 x 10 30 cm 3 ( . (6) 



m DM / 

Neglecting evaporation for the moment, the present WIMP annihilation rate is given by 

r = l -A®N 2 = l -C & tanh 2 (t e /r E ) , (7) 

where t ~ 4.5 billion years is the age of the solar system, and te = (C & A & ) is the time 
required to reach equilibrium. The annihilation rate is maximized when it reaches equilibrium 
with the capture rate. This occurs when 

t Q /r E > 1 . (8) 
In the case in question, this condition will be met so long as 

, n /lGeV\ 1/2 /10- 40 cm 2 \ 
(av) > 3 x 10" 30 cm 3 /s . (9) 

V m DM / V °H / 

Once equilibrium is reached, the final annihilation rate (and corresponding neutrino flux and 
event rate) has no further dependence on the dark-matter particle's annihilation cross section. 

For WIMPs with masses in the range being considered here, the process of WIMP evapora- 
tion from the Sun could also be potentially important. Ref. [28[ found an approximate timescale 
for WIMP evaporation given by (see also Ref. [291]): 

E & W Kr(K"WGeV) + 4) B -l ( ) . (10) 

V5xl0- 39 cm 2 y v ; 

The WIMP annihilation rate in the presence of evaporation is 29l 
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FIG. 1: The factor by which the WIMP annihilation rate in the Sun is suppressed as a result of WIMP 
evaporation. For each WIMP mass, we used a spin-independent elastic scattering cross section near 
the middle of the DAMA region (see the upper frame of Fig. [3]). For WIMPs heavier than 4 GeV, the 
effect of evaporation is negligible. 
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with a = {1 + (E q t e /2) 2 } ' . For the lightest range of WIMP masses in the DAMA region 
(3-4 GeV), the timescale for the process of WIMP evaporation can be comparable to the time 
required to reach capture-annihilation equilibrium, E°te ^ 1, and thus may potentially reduce 
the annihilation rate. In Fig. [T], we show the suppression of WIMP annihilation resulting from 
evaporation as a function of the WIMP's annihilation cross section, obtained from the ratio of 
Eq. (jlip and Eq. ([7j). For WIMPs heavier than ~ 4 GeV, evaporation plays a negligible role. 
For lighter WIMPs, however, the annihilation rate and corresponding neutrino flux may be 
suppressed, depending on the magnitude of the annihilation cross section. In our calculations, 
we will use a dark-matter annihilation cross section of av = 3 x 10 -26 cm 3 /s, which is generally 
expected for a thermal relic in the absence of resonant annihilation, coannihilations or strong 
s-wave suppression. If the annihilation cross section is smaller than this value, the effects 
of WIMP evaporation on the annihilation rate and corresponding neutrino flux may be more 
pronounced, as shown in Fig. [TJ 



III. LIMITS FROM MUON PRODUCTION IN SUPER-KAMIOKANDE 

As the WIMPs annihilate, they can generate neutrinos through a wide range of channels. 
Annihilations to bottom quarks, charm quarks, and tau leptons each generate neutrinos in 
their subsequent decays. In some models, WIMPs can also annihilate directly to neutrino 
pairs. Once produced, neutrinos travel to the Earth where they can be detected. The spectra 
of muon neutrinos and antineutrinos at the Earth from WIMP annihilations in the Sun is given 
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by 
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where C & is the WIMP capture rate in the Sun, F Eq is the non-equilibrium suppression factor 
(« 1 for capture-annihilation equilibrium), Z^es is the Earth-Sun distance and the bracketed 
quantities are the neutrino and antineutrino spectra from the Sun per annihilating WIMP. 



Due to — v T vacuum oscillations and MSW enhanced v e oscillations [30TI in the Sun, the 



muon neutrino flux is approximately given by the average of the u e and v T components, 
leading to the factor of 1/3 in Eq. (|12p . The oscillations of anti-electron neutrinos are MSW 
suppressed, however, leading to a flux of anti-muon neutrinos which is the average of the 
and v T components |31| . 

Muon neutrinos produce muons in charged-current interactions with nuclei inside of and 
around the detector volume of Super-Kamiokande. In the analysis of the Super-Kamiokande 



collaboration [22J, upward-going muon tracks extending 7 meters or more within the inner 
detector were counted. Such events are produced through dark-matter annihilations in the Sun 
at a rate given by 
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Here, a v {E Uy ) » 8 x 10~ 39 cm 2 x [E u /{GeV)} and a D {E^) « 3 x 1(T 39 cm 2 x [E p /{GeV)} are the 
neutrino-nucleon and antineutrino-nucleon charged-current interaction cross sections, (1 — y) is 
the fraction of the neutrino's (or antineutrino's) energy which goes into the muon, v4n e t is the 
effective area of the detector, Na is Avogadro's number, which for water is the number density 
of protons, and R^Ey) is the range a muon travels before losing its energy. For the dimensions 
of the inner detector, we take A^et — 900 m 2 and a height of 36.2 meters, for a total target 
volume of 32,500 m 3 , or 32.5 metric tons. To account for the 7-meter cut applied in the Super- 
Kamiokande analysis, we substitute the physical muon range (i? M ~ 5 meters xE^/GeV), with 
zero if R^ < 7 meters or otherwise with R^ + (36.2 — 7) meters. The factor of 1/2 accounts for 
the Sun being below the Super-Kamiokande detector approximately 50% of the time. 

For the injection spectrum of neutrinos produced in WIMP annihilations, we consider the 
following annihilation channels: uu, t + t~, cc, and bb. We neglect annihilation products such 
as muon pairs, light mesons, etc. as they are expected to lose the vast majority of their energy 
before decaying and thus do not produce energetic neutrinos. For the case of annihilations to 
neutrino-antineutrino pairs, we have assumed that equal quantities of each flavor are produced. 
If instead all of the neutrinos and antineutrinos produced were of electron (muon or tau) flavor, 
the rate would be 27% smaller (14% larger). 

For annihilations to tau pairs, we include the semi-leptonic decays r — > fiuu, euu, as well 
as from the hadronic decays r — > nu, Ku, mru, and t\i\i\v. For charm and bottom quarks, 
only semileptonic decays contribute. We model the neutrino energy spectra coming from these 
decays using updated versions of the formulae in Ref. 32| including tau and electron neutrinos. 
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For the hadronic r decays we use approximate formulae that reproduce reasonably well the 
results obtaining using HERWIG [33J. We neglect the possible production of quarkonia near 
the bb threshold and instead use the results for open-flavor semi-leptonic decay; as the neutrinos 
from 6-decay do not impose strong constraints we believe this approximation is sufficient. 

In Fig. [2J we plot the rate of events at Super-Kamiokande predicted from annihilating 
WIMPs. In the upper frame, we consider spin-independent elastic scattering with nuclei, which 
in Eq. (pQ) corresponds to o"h = a p,si an d cr He = 16 a Pt si- In the lower frame, we consider spin- 
dependent scattering with oh = cr p ,SD and o"H e = (we neglect the tiny contribution from 
spin-dependent scattering from 3 He). For comparison, we also show as thin (blue) lines in 
these frames the rate neglecting the effects of WIMP evaporation. 

We can translate our predicted rate in Super-Kamiokande to a limit on the WIMP elastic 
scattering cross section as a function of mass and the dominant annihilation channel. According 
to the analysis of the Super-Kamiokande collaboration 22|], 170 events were observed within an 
angular window of radius 36° centered around the Sun, over 1679.6 live days. Comparing this 
to the expected rate from atmospheric neutrinos of 185 events, this leads to a 2a upper limit on 
the contribution from WIMP annihilations of approximately: [170 + 2V170] — 185 ~ 11 events. 
We find a similar result is a somewhat smaller angular window is considered. 

In Figs. [3] and HI we apply this limit and plot our results in the (cr, tb^m) plane for the 
cases of spin-independent and spin-dependent scattering. We show both the two-bin analysis 
procedure of Ref. 0] and the full spectral analysis of the DAMA data both with and without 
the 2 — 2.5 keVee bin at threshold [9]. We include in these plots the constraint coming from 
demanding that the total rate observed by the DAMA collaboration not exceed that predicted 
by a given WIMP mass and cross section. The form factors and couplings required for the 



calculation of the spin-dependent DAMA allowed region were obtained from Refs. [37J, [38 
We find that if the WIMP annihilates to neutrinos or taus a significant fraction of the time, 
existing data from Super-Kamiokande closes a large fraction of the DAMA window. In the 
case of spin-dependent scattering, WIMPs in the DAMA region which annihilate to neutrinos, 
tau leptons, or charm quarks or bottom quarks with any significant probability are ruled out 
by Super-Kamiokande. The only possible exception occurs near wdm ~ 2.6 — 3.1 GeV, where 
evaporation enables the WIMP to evade detection by Super-Kamiokande. 



IV. LIGHT NEUTRALINO DARK MATTER 

Within the context of the R-parity conserving minimal supersymmetric standard model 
(MSSM), a neutralino (if the lightest supersymmetric particle) in the mass range of the DAMA 
window would very likely be overproduced in the early Universe [39]. The main reason for this 
conclusion is that LEP constraints force the MSSM Higgs bosons, charginos and sfermions to 
be heavier than ~ 100 GeV, in which case they are unable to efficiently mediate the process 
of neutralino annihilation or to participate in coannihilations. In supersymmetric models with 
extended Higgs sectors, however, this conclusion can be evaded. In the next-to-ninimal super- 
symmetric standard model (NMSSM), for example, it has been shown that very light neutralinos 
(~1-10 GeV) can be thermally produced in acceptable quantities |!o| (see also Ref. [Il|]). This 
is made possible by the exchange of a relatively light pseudoscalar Higgs boson. Intriguingly, if 
the lightest Higgs scalar is very light and singlet-like, light neutralinos are naturally predicted 
to possess an elastic scattering cross section near or within the DAMA window. 

Being Majorana fermions, neutralino annihilation to fermions is chirality suppressed [42], 
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FIG. 2: The number of upward-going muon events per year in Super-Kamiokande from WIMPs 
annihilating in the Sun as a function of mass for an elastic scattering cross section with protons 
of 10 -40 cm 2 , assuming 100% annihilation to the indicated channel. The upper and lower frames 
correspond to spin-independent and spin-dependent couplings, respectively. In each frame, the thin 
(blue) lines extending to the left denote the results neglecting the effects of WIMP evaporation. See 
text for more details. 



leading to av oc rn? f /m 2 . For neutralinos in the mass range under consideration here, an- 
nihilations proceed almost entirely to combinations of bb, t + t~ and cc. Above the bottom 
quark threshold, neutralino annihilations mediated via Higgs exchange produce bottom quarks 
approximately 3m|/ (3m 2 . + 3m 2 + m 2 ) ~ 90% of the time. For neutralinos near the b mass, 
however, the phase space for annihilations to bb is reduced and the fraction of annihilations 
producing tau leptons and charm quarks is enhanced. Below the bb threshold, at least 40% 
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FIG. 3: The limit on a light WIMP's spin-independent elastic scattering cross section with nuclei from 
Super-Kamiokande for various choices of dominant annihilation modes. The upper frame contains the 
DAMA allowed region as calculated in the two-bin analysis of Ref . [fj , while the lower frame uses the 
full spectral analysis with (dark hatched region) and without (light hatched region) the 2 — 2.5 keVee 



bin. Also shown are the limits from the CDMS [10], CRESST [11], CoGeNT [TJ], and XENON[13] 
collaborations. The dotted green lines are the constraints derived by demanding that the predicted 
total rates predicted by a given WIMP candidate do not exceed those observed by DAMA at the 2a 
level in any energy bin. 



of neutralino annihilations proceed to t + t~ (and a considerably larger fraction if the Higgs' 
couplings to down-type fermions are enhanced by tan/5). 

Comparing this to the results shown in Figs. [3] and HI we can place constraints on light neu- 
tralinos as the source of the DAMA signal. Considering the case of spin-independent scattering, 
neutralinos above the bb threshold (4.2 GeV < m x o < 8 GeV) annihilate mostly to bottom 
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FIG. 4: The limit on a light WIMP's spin-dependent elastic scattering cross section with nuclei from 
Super-Kamiokande for various choices of dominant annihilation modes. The upper frame contains 
the DAMA allowed region as calculated following the two-bin analysis of Ref. [5J|, while the lower 
frame uses the full spectral analysis with (dark hatched region) and without (li ght hatched region) the 
2 — 2.5 keVee bin. Also shown are the limits from the CDMS 



3J, CRESST [U, XENON [35| and 



COUPP 36] collaborations. The dotted green lines are the constraints derived by demanding that the 
predicted total rates predicted by a given WIMP candidate do not exceed those observed by DAMA 
at the 2<T level in any energy bin. 



quarks and are not able to be constrained by Super-Kamiokande beyond the corresponding 
limits from XENON, CDMS, and CoGeNT. Below the bb threshold (3 GeV < m x o < 4.2 GeV), 
however, many or most neutralino annihilations produce tau pairs, allowing Super-Kamiokande 
to exclude this range of neutralino masses. For spin-dependent scattering, the constraints from 
Super-Kamiokande are even more severe. It is difficult to imagine any neutralino that could 
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generate the observed DAMA signal through spin-dependent scattering without being in conflict 
with the constraints we have presented in this paper (the possible exception being a neutralino 
with a ~ 3 GeV mass, which could efficiently evaporate in the Sun, thus suppressing the 
annihilation rate). 

V. CONCLUSIONS 

In this paper, we have discussed the indirect detection of light WIMPs in the several 
GeV range by their capture in the Sun and subsequent annihilation to neutrinos. Such 
dark-matter particles can potentially accommodate the annual modulation signal observed by 
DAMA/LIBRA and the null results of other direct-detection experiments, once the effects of 
channeling are accounted for. We consider the constraints that can be placed by the Super- 
Kamiokande experiment on WIMPs in the DAMA window for the cases of prompt dark-matter 
annihilation to neutrinos, and secondary neutrinos produced by decays of tau leptons and charm 
and bottom quarks. The constraints found are significant, and impose stringent constraints on 
the DAMA allowed parameter space, particularly in the low mass end of the DAMA region. 
For the case of spin-independent elastic scattering of WIMPs with nuclei, dark-matter parti- 
cles that annihilate to tau leptons or neutrinos a significant fraction of the time are excluded 
by Super-Kamiokande measurements. For spin-dependent scattering, the constraints are more 
severe; any annihilation fraction to neutrinos, tau leptons, or charm or bottom quarks above 
the 10~ 2 level is ruled out. 
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